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Topology and chiral symmetry in finite temperature QCD* 

Robert G. Edwards a , Urs M. Heller a , Joe Kiskis b and Rajamani Narayanan 
a SCRI, The Florida State University, Tallahassee, FL 32306-4130, USA 
b Dept. of Physics, University of California, Davis, CA 95616, USA 
c American Physical Society, One Research Road, Ridge, NY 11961, USA 

We investigate the realization of chiral symmetry in the vicinity of the deconfinement transition in quenched 
QCD using overlap fermions. Via the index theorem obeyed by the overlap fermions, we gain insight into the 
behavior of topology at finite temperature. We find small eigenvalues, clearly separated from the bulk of the 
eigenvalues, and study the properties of their distribution. We compare the distribution with a model of a dilute 
gas of instantons and anti-instantons and find good agreement. 
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1. Some properties of the overlap Dirac 
operator 

The overlap Dirac operator is constructed out 
of a Wilson-like Dirac operator (the usual Wilson 
Dirac operator is used for this talk) with a Jarge 
negative mass |l],|j : 



D(H) = - [1 + M + (1 - M)7 5 e(# w (m))] 



(1) 



Here H w (m) — 7s.Dwilson( — m). This form 
insures that topology is seen by the overlap 
fermions, since Q = tre(H w )/2. The quark mass 
is proportional to /i (— 1 < fx < 1) for small \x. 
The propagator for external fermions is given 

by 



(2) 



i.e. it has a contact term subtracted. The prop- 
agator is then chiral in the massless case. 

In many cases, it is more convenient to use the 
hermitian version H (fj,) — j 5 D(fi). The massless 
version satisfies, 



{if o (0), 75 } = 2ff o 2 (0). 



(3) 



It follows that [i/^(0),75] = 0, i.e. the eigen- 
vectors of -ffo(O) can be chosen as chiral. Since 



iT o 2 ( M ) = (l-/i 2 )# o 2 (0) + 



/' 



(4) 



this holds also for the massive case. 

Each eigenvalue < A 2 < 1 of -ff 2 (0) is doubly 
degenerate with opposite chirality eigenvectors. 
In this basis, H (n) and D(u) are 2x2 block 
diagonal. In D{jj) the blocks are 



(1 - M )A 2 + M 



-(i-aOAVI 

where 





A 2 



(i - /i)Av / r 



- A 2 



75 = 



-1 



(5) 



(6) 



For a gauge field with topological charge Q ^ 
0, there are, in addition, \Q\ exact zero modes 
with chirality sign(Q) paired with eigenvectors of 
opposite chirality and eigenvalue 1. These are 
also eigenvectors of H (u) and D(fi): 



(7) 



^(m)zcio sector 

depending on the sign of Q 
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This is the Ginsparg- Wilson relation, which is of- 
ten used these days to prove good chiral proper- 
ties of the overlap Dirac operator. 



* Presented by U. M. Heller at Lattice '99. 



2. Small eigenvalues and the chiral con- 
densate 

In the chiral eigenbasis of i? 2 (0), the external 
propagator takes the block diagonal form with 
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2x2 blocks 



ffi( l-X 2 ) -A^/l - A 2 
^A^I^A 2 ^(l-A 2 ) 



(8) 



In topologically non-trivial background fields, 
there are \Q\ additional blocks 



f 





f 



(9) 



depending on the sign of Q. 

We thus find in a fixed gauge field background, 



{^){{U}) 



191 + 1 v 

fjV v ^ 



V ^ A 2 (l 



2/i(l — A 2 ) 
M 2 )+M 2 



(10) 



and averaged over gauge fields, we get the con- 
densate. It is dominated by the small (non- 
zero) eigenvalues, and in the thermodynamic 
limit where the first term vanishes, it is given by 
the density of eigenvalues at zero p(0 + ). 

The close connection between the distribution 
of small eigenvalues and the condensate moti- 
vated us to study the small eigenvalues of the 
overlap Dirac operator in quenched QCD as the 
temperature is raised above the deconfining tran- 
sition temperature T c . A similar study for T = 
can be found in j^j. 

3. Small eigenvalue distribution in 
quenched QCD above T c 

We focus on the non-zero eigenvalues and dis- 
play their distributions for SU(2) and SU(3) 
at various temperatures and spatial volumes in 
Figs. [l]-[3| All figures show a region around 
A = 0.05 with very few eigenvalues, or even none. 
There is also a concentration of small eigenval- 
ues below this value except in the cases with the 
smallest volumes and the largest f3's. Fig. || shows 
that these very small non-zero 

In the remainder, we shall focus on the small 
modes A < 0.05. We denote the number of zero 
and small nonzero eigenvalues with chirality ± 
by n± and the total number by n = n + + n_. 
The topological charge is given by Q = n + — n_. 
Our results for the various ensembles are given in 
Tables and |. We see that for fixed /3, (n) /V 



SU(2) (3=3.3 B d x4 
T « T„ 




SU(3) (3=3.4 B d x4 
T = 1.4 1, 



SU(2) (3=3.4 16 a x4 
T « 1.4 T 



SU(3) (3=3.5 16 a x4 
T = 3 T 
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Figure 1. Small eigenvalue distribution for var- 
ious SU(2) ensembles for temperatures above T c 
on lattices with N t — 4. 




SU(3) (5=5.75 16 3 X4 
T = 1.13 T r 



SU(3) /S=5.75 8 J X4 
T = 1.13 T„ 



SU(3) p = 5.B5 16 X4 
T = 1.38 T„ 



SU(3) = 5.85 8 J X4 
T = 1.38 T„ 



0.0 0.1 0.8 0.0 0.1 0.2 

X X 

Figure 2. Same as Fig. |l| but for SU(3). 
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Figure 3. Small eigenvalue distribution in various 
topological charge sectors for SU(3) with [3 = 5.75 
on a 16 3 x 4 lattice. 

and (Q 2 ) /V seem to remain finite and nonzero in 
the large volume limit, but they drop quickly as 
/?, and hence the temperature, is increased. 



4. Modeling by a dilute instanton 
instanton gas 



anti- 



Looking in more detailed at the small modes, 
we find that their number n is roughly Poisson 
distributed P(n, (n)) = {ri) n e~^ /n\, and in par- 
ticular, the average and variance are approxi- 
mately equal. Also for fixed n, n + and n_ are 
roughly binomially distributed. These observa- 
tions are consistent with interpreting the small 
modes to be due to a dilute gas of instantons and 
anti-instantons with n+ and n_ their numbers. 
n — \Q\ of the would-be zero modes overlap and 
mix to give small eigenvalues, while \Q\ exact zero 
modes remain. 

At finite temperature, instantons fall off expo- 
nentially, and so do the fermionic zero modes as- 
sociated with them. We consider a toy model 
of randomly (Poisson and binomially) distributed 
instantons and anti-instantons, inducing interac- 
tions of the form /ioe _d ( J ' J ^ D between the would- 
be zero modes of every instanton - anti-instanton 
pair with separation d(i,j). Like sign pairs 
are assumed to have no interactions. This toy 



Table 1 

(n) and (Q 2 ) per 8 3 spatial volume and (n)/a n 
with a n the variance of n, for our SU(2) ensem- 
bles. 



volume 




8 a x 4 
2.3 2.4 


16 a x 4 
2.4 2.5 


(n)/V 

(Q 2 )/v 

(n)/a n 


1.7 0.29 
1.7 0.31 
0.98 1.09 


0.25 0.054 
0.25 0.053 
0.93 0.99 



Table 2 

Same as Table | but for SU(3). 



volume 

P 


8 3 x 4 
5.75 5.85 


12 3 x 4 
5.75 


16 3 x 4 
5.71 5.75 5.85 


(n)/V 

(Q 2 )/v 

(n)/a n 


0.32 0.058 
0.31 0.068 
1.05 0.90 


0.28 
0.28 
0.92 


0.63 0.31 0.051 
0.64 0.33 0.049 
1.15 1.02 0.83 



model reproduces the qualitative features of the 
small eigenvalue distributions. Preliminary esti- 
mates give D w 1/(2T) 0. 

In conclusion, we find that in quenched QCD 
above the deconfining transition temperature 
topology, as manifested by exact zero modes, per- 
sists. Furthermore, a finite density of small eigen- 
values, separated from the bulk of the eigenval- 
ues, remains. The properties of their distribu- 
tion is well modeled by a dilute gas of instantons 
and anti-instantons. Our observation might be 
a quenched artefact, since in full QCD the finite 
temperature transition is believed to be driven by 
chiral symmetry restoration. 
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